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Abstract. Neutrino Deep Inelastic Scattering (DIS) on nuclei is an essential process to constrain 
the strange quark parton distribution functions (PDF) in the proton. The critical component on the 
way to using the neutrino DIS data in a proton PDF analysis is understanding the nuclear effects in 
parton distribution functions. We parametrize these effects by nuclear parton distribution functions 

^ ' (NPDF) and we use this framework to analyze the consistency of neutrino DIS data with other 

— i ' nuclear data. 
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Ph; introduction 

^ . Any prediction for a process measured at a hadron collider such as the LHC, involves 

^^ ' necessarily parton distribution functions (PDFs). Because of how crucial the knowledge 

of PDFs is, there are many groups that perform and update global analyses of PDFs 
^ . protons [1, 2, 3] and nuclei [4, 5]. Although not often emphasized, nuclear effects are 

CN ! present also in the proton PDFs analysis as a number of experimental data is taken on 

^D • nuclear targets. Mostly though, the nuclear targets used in the proton analysis, are made 

^T) ■ of light nuclei where nuclear effects are generally small. An important exception is the 

r^ ! neutrino DIS data which is taken on heavy nuclei such as iron or lead and is sensitive to 

O I the strange quark content of the proton. A knowledge of the strange quark PDF has an 

influence on precise measurements at the LHC such as W- or Z-boson production. 

In order to include the neutrino DIS data in a global fit to determine proton PDF, we 
have to apply a nuclear correction factor. The nuclear correction factor can be obtained 
^ ■ either from a specific model of nuclear interactions or from an analysis of nuclear parton 

^ . distribution functions (NPDF) based on experimental data. 

Here, we present a framework for a global analysis of nuclear PDFs at next-to-leading 
order in QCD closely related to the CTEQ framework for proton PDFs. We analyze and 
compare the nuclear correction factor obtained from the usual charged lepton DIS and 
Drell-Yan (DY) data to the one from the neutrino DIS data mainly from the NuTeV 
experiment. 

NUCLEAR PDF 

To determine parton distribution functions from experimental data, we parametrize the 
jc-dependence of PDFs at an input scale Qo and use the DGLAP equations to change 
the scale of experimental data to Qq in order to perform a fit to the data. The global 
NPDF framework, we use to analyze charged lepton DIS and DY data and neutrino DIS 



data, was introduced in [6]. The parameterizations of the nuclear parton distributions of 
partons in bound protons at the input scale of Qq = 1.3GeV 

xfkix^Qo) = co^^-'(l -x)^-2e^-3-(l +e''xy^ , (1) 

where k — Uv,dvjgjU + J, 5, s and 

d{x, Qq) juix, Qq) = CQX'^ {l-x)'^ + {\+ C3x) {l-x)'\ (2) 

are a generalization of the parton parameterizations in free protons used in the CTEQ 
proton analysis [7]. To account for different nuclear targets, the coefficients Ck are made 
to be functions of the nucleon number A 

c,,^Ck{A) = Ck,o + Ck,i{l-A-'^-^), A:={1,...,5}. (3) 

The proton PDF in this framework are obtained as a limit A — )> 1 and are held fixed at 
values obtained in the analysis [7]. From the input distributions, we can construct the 
PDFs for a general (A, Z) -nucleus 

t^'\x.Q) = J f!'\x,Q) + ^^^ff\x,Q), (4) 

where we relate the distributions of a bound neutron, /"' {x, Q) , to those of a proton by 
isospin symmetry. 

In the analysis, the same standard kinematic cuts Q > 2GeV and W > 3.5GeV were 
applied as in [7] and we obtain a fit with ;f ^/dof of 0.946 to 708 data points with 32 free 
parameters (for further details see [6]). 

The nuclear effects extracted in the form of NPDF are usually presented in the form 
of nuclear correction factors. We discuss two nuclear correction factors in the following 
where both are related either to the DIS structure function F2 in the charged-current (CC) 
vA process 

or to the DIS structure function F2 in the neutral-current (NC) l^A process 

The superscript '0' stands for using the free nucleon PDFs /f '"(jc, Q) in Eq. (4). 

In Fig. 1 (solid line), we show how the result of our global analysis of charged 
lepton data translates into these nuclear correction factors and how it compares with 
experimental data. As first observed in [8], i\\QR^(^{F2\x, Q^) correction factor calculated 
using Eq. 5 with parton densities from the fit to the charged lepton nuclear data, does not 
describe the NuTeV data well which raises the question if including neutrino DIS data in 
the global analysis corrects this behavior without spoiling the R^^^{F2;x, Q^) correction 
factor which fits the charged lepton DIS and DY data well. 



TABLE 1. Summary table of a 


family of com- 


promise fits. 
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FIGURE 1. Nuclear correction factors R^^{F2;x,Q ) (left) and R(j^{F2;x,Q ) (right) for compromise 
fits with different weights of the neutrino DIS data at the scale Q^ — 5GeV^. 



NEUTRINO DIS 

In order to analyze the apparent discrepancy between the nuclear correction factor 
R'^^{F2\x, Q^) from the fit to charged lepton data and the neutrino charged current DIS 
data, we have set up a global analysis where we used exclusively the neutrino DIS 
cross-section data coming from NuTeV and Chorus experiments taken on iron and lead 
respectively. Here we applied the same kinematic cuts as in the first analysis of the 
charged lepton data and we obtain a fit to 3134 neutrino DIS cross-section data points 
with x^ /dof of 1.33 with 34 free parameters (for further details see [9]). A global fit 
to neutrino DIS data describes the data for the charged current nuclear correction factor 
R'^(j{F2',x, Q^) well and does a poor job to describe the neutral current correction factor 
especially at low and intermediate Bjorken x. In order to find nuclear correction factors 
in agreement with both charged lepton and neutrino data, we perform a combined global 
analysis of both data sets. We introduce an additional parameter in our analysis, the 
weight of the neutrino data set w, which should prevent the neutrino DIS data set to 
dominate the global fit based only on the number of data points taken. The weight w 
enters the calculation of the combined ;t^ as 



X 



/±A data 



zf + 



vA data 



(7) 



and it can be used to interpolate between the two different global fits (w = results in the 
fit to charged lepton data only and w = oo stands symbolically for the fit only to neutrino 
data). Varying the weight w, we try to find a compromise fit which would describe both 



charged lepton and neutrino data well. We list the resulting j} for the compromise fits 
with weights w = 0, 1/7, 1/2, l,oo in Tab. 1 and we show the nuclear correction factors 
B^^{F2\x^ Q^) and R^(j{F2;x, Q^) for the same family of compromise fits in Fig. 1. 

In order to decide on how well the compromise fits describe the data we use the X^ 
goodness-of-fit criterion used in [10, 2]. We consider a fit a good compromise if its X^ 
for both data subsets, the charged lepton DIS and DY data and the neutrino DIS data, is 
within 90% confidence level of the fits to only charged lepton or neutrino data. 

We define the 90% percentile 1^90 used to define the 90% confidence level, by 

P{x\N)dx^ = 0.90, (8) 



, (y2\N/2-lg-X^/2 

where A'^ is the number of degrees of freedom and P{x ,N) = JN/2r(Aj/'y\ — i^ the 

probability distribution. We can assign a 90% confidence level error band to the x^ 
of the fits to the charged lepton DIS and DY data and to the neutrino DIS data 

;tf±^ = 638 + 45.6, Xva = '^^92+\3S. (9) 

Comparing the results of the compromise fits with different weights, listed in Tab. 1, we 
see that none of the compromise fits are compatible with both 90% confidence level lim- 
its given in Eq.9. As detailed in [9], not even relaxing the condition to compare against 
the 99% confidence level limit helps to finding a suitable compromise fit. Moreover, we 
show in [9] that the effect is related to the precise neutrino DIS data from NuTeV. 



CONCLUSION 

After performing a thorough global NPDF analysis of the combined charged lepton and 
neutrino data, we find that there is no good compromise description of both the data 
sets simultaneously. The differences are most pronounced in the low and intermediate x 
regions where the neutrino DIS (NuTeV) do not show a strong shadowing effect as the 
charged lepton data do. The inability to describe all data by one consistent framework 
indicates the existence of non-universal nuclear effects or unexpectedly large higher- 
twist effects. 
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